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We report on the Dalitz plot analyses of i? — > D+tt tt , 5+ 
plot method and the most recent BABAR. results are discussed. 



and D, 



The Dalitz 



1. Introduction 

Dalitz plot analysis is an excellent way to study 
the dynamics of three-body decays. The decays under 
study are expected to proceed predominantly through 
intermediate quasi-two-body modes [ij and experi- 
mentally this is the observed pattern. Dalitz plot 
analysis can explore resonances properties, measure 
corresponding magnitudes and phases, and investigate 
corresponding CP asymmetries. 

We report on the Dalitz plot analyses of 
D"'"7r~7r~, — > 7r''"7r~7r'*" and Df — > 7r''"7r~7r'*'. 
The data were collected with the BaBAR detector at 
the PEP-II asymmetric-energy e~^e~ storage rings at 
SLAC. A detailed description of the BaBAR detector 
is given in Ref. 01 . Monte Carlo (MC) simulation is 
used to study the detector response, its acceptance, 
background, and to validate the analysis. 



rj 2. Dalitz Plot Analysis 



In the decay of a spin-0 particle ^ to a final state 
composed of three pseudo-scalar particles (abc), the 
differential decay rate is generally given in terms of 
the Lorentz-invariant matrix element M by 



dxdy 



\M\ 



(1) 



where is the A particle mass, x and y are the 
invariant masses-squared of the ab and ac pairs, re- 
spectively. The Dalitz plot provides a graphical rep- 
resentation of the variation of the square of the matrix 
element, |A4p, over the kinematically accessible phase 
space {x,y) of the process. 

The distribution of candidate events in the Dalitz 
plot is described in terms of a probability density func- 
tion (PDF). The PDF is the sum of signal and back- 
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ground components and has the form: 

Bix,y) 



PDF(a;,y) = /bg 



+ (l-/bg) 



/dp ^(^^ y)d'Xdy 

[S{x,y)®n]e{x,y) 



/dp [Six,y) iE''R]e{x,y)dxdy' 

(2) 



where 



• The integral is performed over the whole Dalitz 
plot (DP). 

• S {— |A1p) and B describe signal and back- 
ground amplitude contributions, respectively. 

• The S{x, y) ® TZ is the signal term convolved 
with the signal resolution function, the signal 
resolution can be safely neglected for broad res- 
onances. 

• /bg is the fraction of background events. 

• e is the reconstruction efficiency. 

An unbinned maximum likelihood fit to the Dalitz 
plot is performed in order to maximize the value of 



C= II PDF(a;„y, 



(3) 



z=l 



with respect to the parameters used to describe S, 
where Xi and yi are the values of x and y for event i 
respectively, and A*"event is the number of events in the 
Dalitz plot. In practice, the negative-log-likelihood 
(NLL) value 



NLL 



-InC 



(4) 



is minimized in the fit. 

The isobar model formulation is used, in which the 
signal decays are described by a coherent sum of a 
number of two-body amplitudes. The total decay ma- 
trix element A4 is given by 



M 



^Pke 

k 



(5) 



2 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



where pk and ipk are the magnitudes and phases of the 
fc"^ decay mode respectively, the Ak distributions for 
resonant contribution describe the dynamics of the de- 
cay amphtudes and are written as the product of a rel- 
ativistic Breit-Wigner function, two Blatt-Weisskopf 
barrier form factors [1] and an angular function Q- 

The efficiency-corrected fraction due to the resonant 
or non-resonant contribution k is defined as follows: 



\Pk\ 



' /j3p \ Ak\'^dxdy 



/dp I P '^^^y 



(6) 



The fk values do not necessarily add to 1 because 
of interference effects. The uncertainty on each fk is 
evaluated by propagating the full covariance matrix 
obtained from the fit. 



3. Results from Bj^ar 



3.1. B- 



TT TT 



ISA 



We reconstruct the decays B D^tt tt [6] with 
_D+ K^TT^TT^ based on a sample of about 383 x 
10^ T(4S) BB decays. The resuhing Dalitz plot 
distribution for data is shown in Fig. [TJ The number 
of signal events in the Dalitz plot is 3414 ± 85. The 
background fraction is about 30.4%. The background 
is parametrized using the MC simulation. 



and m^;„(I?7r) are the heavy and light invariant 
masses-squared of the Dtt systems, respectively. The 
distributions in Fig. [2] show good agreement between 
the data and the fit. 

We have performed some tests with the broad reso- 
nance Dq^ excluded or with the of the broad res- 
onance replaced by other quantum numbers. In all 
cases, the NLL values are significantly worse than that 
of the nominal fit. We conclude that a broad spin-0 
state Z)q° is required in the fit to the data. 

The total branching fraction of the B~ D+tt^tt^ 
decay is measured to be 



B{B- 



) ^ (1.08 ±0.03 ±0.05) X 10' 



where the first error is statistical and the second is sys- 
tematic. The systematic error on the measurement of 
the total B- D+ TT TT branching fraction is due 
to the uncertainties on the number of B'^B~ events, 
the charged track reconstruction and identification ef- 
ficiencies, the —t K^TT^TT^ branching fraction, the 
background shape and the fit models. 

The mass and width of I?2° are determined to be: 



(2460.4 ± 1.2 ± 1.2 ± 1.9) MeV/c^ and 
(41.8 ± 2.5 ± 2.1 ± 2.0) MeV, 



respectively, while for the Dq'^ they are: 
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Figure 1: Dalitz plot for B 



r)*o 



r)*0 



The nominal fit model contains 
Z?*(2007)°, B* and "P-wave non-resonant contribu- 
tions. The 5- wave and 2?- wave non-resonant contri- 
butions are found to be small and can be ignored. 
Replacing the D*(2007)° by a Dn 5- wave [7| state is 
found to give a much worse fit. 

Fig. Hi, [2d andU; show the ml^.^^{Dn), TO^^^(D7r) 
and m^(7r7r) projections respectively. Here m^-^,^^{DTT) 



mo-o = (2297 ± 8 ± 5 ± 19) MeV/c^ and 
r^,o = (273 ± 12 ± 17 ± 45) MeV, 

where the first and second errors reflect the statistical 
and systematic uncertainties, respectively, the third 
one is the uncertainty related to the fit models and 
the Blatt-Weisskopf barrier factors. 

We have also obtained exclusive branching fractions 
for and production: 

B{B- Dfir-) X B{Df -> ZJ+tt") 

= (3.5 ± 0.2 ± 0.2 ± 0.4) X 10"'' and 

B{B- D*°Tr-) X B{D*° D+n-) 

= (6.8 ± 0.3 ± 0.4 ± 2.0) X 10"*. 

The relative phase of the scalar and tensor ampli- 
tude is measured to be 



-2.07 ±0.06 ±0.09 ±0.18 rad. 



The systematic uncertainties that affect the results 
of Dalitz plot analysis are due to the background 
parametrization, the selection criteria, and the pos- 
sible fit bias. 

Our results for the masses, widths and branching 
fractions are consistent with but more precise than 
previous measurements performed by Belle Q. 
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Figure 2: Result of the nominal fit to the data: projec- 
tions on (a) m^^^{Dn), (b) rri^^^{D'jT) and (c) m^(7r7r). 
The points with error bars are data, the solid curves rep- 
resent the nominal fit. The shaded areas show the 
contribution, the dashed curves show the Dq" signal, the 
dash-dotted curves show the _D*(2007)'' and B* signals, 
and the dotted curves show the background. 



3.2. B+ 



TV TT TT 



We fit the B" and B'^ samples independently in or- 
der to determine the CP asymmetry. The total signal 
amplitudes for B~^ and B~ decays are given by 



E 

j 



Cj Aj (^rrij^-^^^, nijjjjjj) , 



(7) 



where rn^^^ and rn^-^^ are the heavy and light invari- 
ant masses-squared of the tt^tt^ systems, respectively. 
The complex coefficients Cj and Cj for a given decay 
mode j contain all the weak phase dependence. Since 
the Aj terms contain only strong dynamics, Aj = Aj. 
The CP asymmetry for each contributing resonance is 
determined from the fitted parameters 



Acp, i — 



(8) 



7r+ based 



We reconstruct the decays i?+ tt^ 
on a sample of about 465 x 10^ T(4S) BB decays. 
We reject background from two-body decays of 
meson and charmonium states {J/ip and ■0(25')) by 
excluding invariant masses (in units of GeV/c^) in the 
ranges: 1.660 < m^+^- < 1.920, 3.051 < m^+^- < 
3.222, and 3.660 < m^+^- < 3.820. The background- 
subtracted Dalitz plot distribution for data is shown 
in Fig. [3l We fit 4335 B candidates to calculate the 
fit fractions and CP asymmetries. The background is 
parametrized using the MC simulation. 
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Figure 3: Background-subtracted Dahtz plot for B 
tt'^tv'tt'^ . The plot shows bins with greater than zero en- 
tries. The depleted bands are the charm and charmonia 
exclusion regions. 

The nominal signal Dalitz plot model comprises 
a momentum-dependent non-resonant component 
and four intermediate resonance states: p°(770)7r*, 
p°(1450)7r±, /2(1270)7r±, and /o(1370)7r±. We do 
not find any significant contributions from /o(980)7r^, 
Xcoi"^; or Xc27r^. The absence of the charmonium 
contributions precludes the extraction of the unitar- 
ity triangle angle 7. We find no evidence for direct CP 
violation. The inclusive CP asymmetry is determined 
to be (-h3.2 ± 4.4 ± 3.1 ^2:0) "^O' where the uncertain- 
ties are statistical, systematic, and model-dependent, 
respectively. 

Projections of the data, with the fit result overlaid, 
as TT^TT^ invariant-mass distributions can be seen in 



4 
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Figure 4: Dipion invariant mass projections: (upper) in 
the p^(J7Q) region; and (lower) in the regions of XcO and 
Xc2- The data are the points with statistical error bars, 
the dark-shaded (red) histogram is the qq component, the 
light-shaded (green) histogram is the BB background con- 
tribution, while the upper (blue) histogram shows the total 
fit result. The dip near 0.5 GeV/c^ in the upper plot is due 
to the rejection of events containing K*^ candidates. 



parametrization, the charged track reconstruction and 
identification efficiencies, and the possible fit bias. 

Our results will be useful to reduce model uncer- 
tainties in the extraction of the CKM angle a from 
time-dependent Dalitz plot analysis of Tr+Tr^Tr". 

3.3. Dt n+iT-n+ (TS 

We reconstruct Z)+ — > 7r+7r~7r+ based on 384 fb~^ 
data sample. We use the tag of £)*(2112)+ Df-/ 
to reduce the combinatorial background. The signal 
PDFs are build using the D+ K+K-Tr+. The re- 
constructed sample contains 13179 events with a 
purity of 80%. The background shape is obtained by 
fitting the Df sidebands. The Dalitz plot distribution 
for data is shown in Fig. [5l 
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FiglH The agreement between the fit results and the 
data is generally good. 

The Dalitz plot is dominated by the /9*'(770) res- 
onance and a non-resonant contribution. The mass 
and width of the /o(1370) are determined to be 
™/o(i370) = 1400 ±40MeV/c2 and rj„(i370) = 300 ± 
80MeV (statistical uncertainties only). We mea- 
sure the branching fractions B{B^ tt^tt^tt^) ~ 
(15.2±0.6±1.2±0.4) X 10-6, B(S± ^ p°{770)n'^) = 
(8.1 ± 0.7 ± l.2t"i) X 10-^ S(B± ^ /2(1270)7r±) = 
(1.57 ± 0.42 ± 0.16+°H) X 10-^ and B(B± ^ 
7r=^7r=^7r=F non - resonant) = (5.3 ± 0.7 ± 0.6^^:5) x 
10~6, where the uncertainties are statistical, sys- 
tematic, and model-dependent, respectively. We 
have made the first observation of the decay 
/2(1270)7r='= with a statistical significance of 6.1a. 

The systematic uncertainties that affect the re- 
sults of Dalitz plot analysis are due to the the num- 
ber of B^B^ events, signal and background PDF 



Figure 5: Dalitz plot for D'^ —* tt^tt tt^. 

For the tt+tt- 5-wave amplitude, we use a model- 
independent partial wave analysis: instead of includ- 
ing the iS-wave amplitude as a superposition of rela- 
tivistic Breit-Wigner functions, we divide the tt+tt" 
mass spectrum into 29 slices and we parametrize the 
iS-wave by an interpolation between the 30 endpoints 
in the complex plane: 

^5-wavc("iTr7r) = Interp(cfc (m^^ )e**'= ("""^ )fc=i, . . ,30 . 

(9) 

The amplitude and phase of each endpoint are free 
parameters. The width of each slice is tuned to get 
approximately the same number of tt+tt- combina- 
tions (~ 13179 X 2/29). In terp olation is implemented 
by a Relaxed Cubic Spline [ll| . The phase is not con- 
strained in a specific range in order to allow the spline 
to be a continuous function. 

The nominal signal Dalitz plot model includes 
/2(1270), /5(770), p(1450) and tttt 5-wave. The Dalitz 
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Figure 6: (a) <S-wave amplitude extracted from the best fit, (b) corresponding <S-wave phase. 
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Figure 7: Dalitz plot projections (points with error bars) and fit results (solid histogram), (a) m^(7r^7r )iow, (b) 
m^(7r^7r~)high, (c) total m^(7r+7r~), (d) mF'{'K^-K^). The hatched histograms show the background distribution. 



plot is dominated by the Dg —>■ (7r"'"7r~)5_wavo7r^ 
contribution. The iS-wave shows, in both amphtude 
and phase, the expected behavior for the /o(980) res- 
onance, and further activity in the regions of the 
/o(1370) and /o(1500) resonances. The iS-wave is 
small in the /o(600) region, indicating that this reso- 
nance has a small coupling to ss. 

The resulting iS-wave tt+tt^ amplitude and phase is 
shown in Fig. [5^ andO). Our results on the amplitude 
and phase of the tt+tt^ iS-wave agree better (within 
uncertainties) with the results from FOCUS than 
those from E791 [l^. The Dalitz plot projections to- 
gether with the fit results are shown in Fig. [71 



Since D+ 



'7r+ and Dt — > K~^K tt"*" have 



similar topologies, the ratio of branching fractions is 
expected to have a reduced systematic uncertainty. 
We therefore select events from the two Dg decay 
modes using similar selection criteria for the D*~^ se- 
lection. The ratio of branching fractions is evaluated 



B{D+ K+K-1T+) 



0.199 ±0.004 ±0.009, 



where the first error is statistical and the second is 
systematic. The systematic uncertainty comes from 
the MC statistics and from the selection criteria used. 
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Conclusion the TT+TT <S-wave is extracted in a model- 

independent way for the first time. 



4. 

We reviewed the recent results from the Dalitz plot 
analyses of B~ — *■ D+7r~7r~, -B+ — *■ 7r+7r~7r+ and 
7r"'"7r~7r"'~ performed by the BaBAR, Collabora- 
tion. The results can be summarized as follows: 

• B- £>+7r-7r- 
We find the total branching fraction of the three- 
body decay: B{B' D+tt-tt-) = (1.08 ± 
0.03 ±0.05) X 10-3. We observe the established 
and confirm the existence of Dq*^ in their 
decays to D~^tt~, we measure the masses and 
widths ofD*° and D*° to be: m^.o = (2460.4 ± 

1.2±1.2±1.9)MeV/c2, T^.o = (41.8±2.5±2.1± 
2.0) MeV, m^.o = (2297±8±5±19) MeV/c^ and 
rjj,o = (273 ± 12 ± 17 ± 45) MeV. 

We measure the branching fractions B{B^ — > 
7r±7r±7r=F) = (15.2 ± 0.6 ± 1.2 ± 0.4) x 
10-6, B{B^ /5"(770)7r±) = (8.1 ± 0.7 ± 
1.2tii) X 10"®> >3iB^ /2(1270)7r±) = 
(1.57 ± 0.42 ± 0.16 1'^il) X 10-^ and B(B± ^ 
TT^TT^TT^ non — resonant) = (5.3 ± 0.7 ± 
0.6l J;5) X 10-6. We observe no significant direct 
CP asymmetries for the above modes, and there 
is no evidence for the decays B^ — > /o(980)7r^, 
XcoTT^, or B^ Xc27r^. 

We perform a high precision measurement of 
the ratio of branching fractions: B{Ds 
Tr+n-TT+)/B{Ds K+K--K+) = 0.199 ± 
0.004 ± 0.009. The amplitude and phase of 
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